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Purpose: A low serum 25-hydroxyvitamin D [25(OH)D] level in the blood has been 
correlated with an increased risk of diabetes mellitus; however, the association be- 
tween serum 25(OH)D level and insulin resistance has not been established in a Ko- 
rean rural population. The aim of this study was to investigate the independent asso- 
ciation between serum 25(OH)D level and insulin resistance in rural Korean adults. 
Materials and Methods: This study used data from the Korean Genome Epidemi- 
ology Study-Kangwha Study. In the 2011 study, 1200 adults completed health ex- 
aminations. In an ancillary study, serum 25(OH)D level was measured in a sub- 
sample (n=813). After excluding those taking vitamin D supplements, a cross- 
sectional analysis was carried out on 807 participants (324 men and 483 women) 
aged 40 to 89 years old. Measured from overnight fasting blood samples, glucose 
and insulin levels were used to calculate the homeostasis model assessment for in- 
sulin resistance (HOMA-IR). Measures of glucose, insulin, and HOMA-IR were 
log-transformed for parametric tests. Results: Serum 25(OH)D level was inversely 
associated with HOMA-IR (P=-0.003, /j=0.039) in a univariate analysis. However, 
the association was not significant after adjustment for sex and age (P=-0.002, 
j?=0.123) or after adjustment for sex, age, body mass index, smoking status, alco- 
hol intake, and regular exercise (P=-0.003, /?=0.247). Conclusion: Our findings 
suggest that vitamin D is not independently associated with insulin resistance in 
Korean men and women. 
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Vitamin D promotes calcium and phosphate absorption in the intestines, maintains 
sufficient concentrations of circulating calcium and phosphate levels, and facilities 
normal mineralization of bone by providing these minerals to bone-forming 
sites. 12 Recently, additional functions for vitamin D have also received greater at- 
tention, especially concerning its effect on insulin resistance. 3 5 Nonetheless, while 
some previous studies reported a significant association between low circulating 
vitamin D levels and insulin resistance, 1 - 6 9 others did not. 10 17 Associations thereof 
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have proven inconsistent and dependent upon the character- 
istics of a particular study population, such as ethnicity, sex, 
and age. 1011 In a prospective cohort study in Denmark, low 
serum 25-hydroxyvitamin D [25(OH)D] status was signifi- 
cantly associated with unfavorable changes in continuous 
markers of glucose homeostasis, but not with incident dia- 
betes. 17 Therefore, it is still unknown as to whether serum 
25(OH)D level has a direct effect on insulin resistance, in- 
dependent of other risk factors. 

The aim of this study was to investigate whether lower 
serum 25(OH)D level is independently associated with in- 
sulin resistance in a rural Korean population. 



MATERIALS AND METHODS 



Study population 

This study used data from the Korean Genome Epidemiol- 
ogy Study (KoGES)-Kangwha study, a rural community- 
based prospective cohort. The KoGES-Kangwha study en- 
rolled 4906 participants aged 40 years or older between 2006 
and 20 1 1 . In the 20 1 1 study, 43 1 new participants completed 
baseline health examinations and 769 participants completed 
follow-up health examinations from July to August. An an- 
cillary biomarker study measured serum 25(OH)D levels for 
a subsample of 8 1 3 participants adopted from the 20 1 1 study. 
We excluded those taking vitamin D supplements (n=4), 
missing serum insulin measurements (h=l), and body weight 
measurements (n=l). Finally, a cross-sectional analysis was 
carried out on 807 participants (324 men and 483 women) 
aged 40 to 89 years old. All participants provided written 
informed consent and the study protocol was approved by 
the Institutional Review Board of Severance Hospital at 
Yonsei University College of Medicine. 

Measurements 

Participants were individually interviewed using standard- 
ized questionnaires to obtain information about their gener- 
al characteristics, medical history, medication use, vitamin 
supplement intake, and lifestyle behaviors, such as smoking, 
alcohol intake, and regular exercise. 18 " 20 Participants were 
classified as either current smokers or current nonsmokers 
and current alcohol drinkers or current nondrinkers. Trained 
interviewers conducted the surveys according to a pre-de- 
veloped protocol and double-checked whether answers were 
missing or inappropriate. 
All participants wore light clothing for anthropometric 



measurements. Standing height was measured up to 0.1 cm 
with a stadiometer (SECA 225, SECA, Hamburg, Germa- 
ny) and body weight was measured up to 0. 1 kg with a dig- 
ital scale (GL-60000-20, Seoul, Korea) according to proto- 
col. Subsequently, body mass index (BMI) was calculated 
as body weight in kg divided by standing height in meters 
squared (kg/m 2 ). Participants were seated for at least 5 min- 
utes before blood pressure measurement, and two measure- 
ments were taken within at least 5 minutes using an automat- 
ic sphygmomanometer (Dinamap 1846 SX/P; GE Health- 
care, Waukesha, WI, USA). If the two measurements differ- 
ed by 10 mm Hg or more for systolic blood pressure (SBP) 
or diastolic blood pressure (DBP), a third measurement was 
taken after 5 minutes, and the last two measurements were 
averaged for analysis. 

Blood samples of all participants were collected from the 
antecubital vein after at least eight hours of fasting. Enzy- 
matic methods were applied to measure total cholesterol, 
high-density lipoprotein cholesterol, and triglycerides (AD- 
VIA 1800; Siemens Healthcare Diagnostics Inc., Deerfield, 
IL, USA). Serum 25(OH)D, currently the best indicator of 
vitamin D store, was measured by radioimmunoassay (DI- 
ASource, Nivelles, Belgium). The intra-assay coefficient of 
variation was 7.3% at 19.1 ng/mL and 8.8% at 10.4 ng/mL. 
The inter-assay coefficient of variation was 7.2% at 6.7 ng/ 
mL and 7.3% at 13.1 ng/mL. Fasting blood glucose level 
was measured using the glucose hexokinase method. Serum 
insulin level was measured in accordance with the radioim- 
munometric method, and the intra-assay coefficient of varia- 
tion was 2.1% at 6.6 uIU/mL and 1.5% at 53.3 uiU/mL. The 
inter-assay coefficient of variation was 6.5% at 14.4 ng/mL 
and 6.1% at 100.4 ng/mL. Additionally, insulin resistance of 
all participants was estimated using the homeostasis model 
assessment for insulin resistance (HOMA-fR), as the prod- 
uct of fasting blood glucose level (mg/dL) and serum insulin 
level (uIU/mL) divided by 405. 21 

Statistical analysis 

We investigated differences between general characteristics 
and concentrations of variables, such as serum 25(OH)D 
with fasting glucose and insulin resistance for men and wom- 
en. Continuous variables were assessed using independent 
t-tests and Wilcoxon rank sum tests, and categorical vari- 
ables were assessed using chi-squared tests. Fasting glucose, 
insulin, and HOMA-IR were skewed to the right, thus log- 
transformed for use in parametric tests. The relationship be- 
tween serum 25(OH)D level and log-transformed HOMA- 
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IR was evaluated using Pearson's correlation coefficients. 
Sex was adjusted when analyzing the total population only. 
Linear regression analyses were used to assess the indepen- 
dent association between serum 25(OH)D and insulin resis- 
tance. First, the analyses were performed unadjusted. Then, 
using multivariate models, we adjusted for 1) sex and age 
and 2) sex, age, BMI, smoking status, alcohol intake, and 
regular exercise. Multiple logistic regression models were 
also used to calculate odds ratios with 95% confidence in- 
tervals of increased insulin resistance (defined as HOMA- 
IR >2.5) at an increment of 10 ng/mL 25(OH)D. 21 SAS ver- 
sion 9.2 (SAS Institute, Cary, NC, USA) was used for all 
analyses, and statistical significance was defined as a two- 
tailed p-vahie <0.05. 



RESULTS 



The general characteristics of all participants are presented 
in Table 1 . This study consisted of 807 participants with a 
mean age of 61.2 years. Mean BMI was 24.5 kg/m 2 , and 
mean SBP and DBP were 114.7 mm Hg and 69.4 mm Hg, 
respectively. Mean serum 25(OH)D level was 25.8 ng/mL 
and mean HOMA-IR was 1.8. The prevalence of diabetes 
mellitus was 15.6%, and current cigarette smoking and al- 



cohol drinking were recorded in 9.2% and 39.0%, respec- 
tively. Additionally, 43.9% of participants reported exercis- 
ing regularly. 

Table 2 presents data on the correlation of serum lipids 
and 25(OH)D and HOMA-IR. Serum 25(OH)D was in- 
versely associated with HOMA-IR in the unadjusted model; 
however, the association disappeared after adjusting for sex 
and age. Supplementary Fig. 1 (only online) presents the lin- 
ear relationship between serum 25(OH)D level and log- 
transformed HOMA-IR separately for men and women. 

Table 3 shows the association between serum 25(OH)D 
level and log-transformed HOMA-IR. Serum 25(OH)D lev- 
el was significantly associated with HOMA-IR (P=-0.003, 
/>=0.039) in the unadjusted model. However, the associa- 
tion was no longer present after adjusting for sex and age 
(P=-0.002, j9=0.123). The association remained non-signifi- 
cant after additional adjustment for BMI, smoking status, 
alcohol intake, and regular exercise (P=-0.002,/>=0.247). 

The association between serum 25(OH)D level and insu- 
lin resistance as a dichotomous measure is shown in Sup- 
plementary Table 1 (only online). Increased insulin resis- 
tance (defined as HOMA-IR >2.5) was not associated with 
serum 25(OH)D. In the unadjusted analysis, the odds of 
having insulin resistance was 0.99 (95% CI, 0.86-1.14) at 
an increment of 10 ng/mL of 25(OH)D. After adjusting for 



Table 1. Characteristics of 324 Men and 483 Women 



Variables 


Total (n=807) 


Men (n=324) 


Women (n=483) 


p value 


Age, yrs 


61.2±9.1 


63.4±8.7 


59.8±9.1 


<0.001 


Height, cm 


158.2±7.9 


165.3±5.4 


153.5±5.2 


<0.001 


Weight, kg 


61.4±9.7 


66.5±9.2 


58.0±8.5 


<0.001 


Body mass index, kg/m 2 


24.5±3.1 


24.3±2.8 


24.6±3.3 


0.168 


Systolic BP, mm Hg 


114.7±15.4 


116.4±14.7 


113.5±15.8 


0.007 


Diastolic BP, mm Hg 


69.4±9.6 


72.4±9.4 


67.3±9.2 


<0.001 


Total cholesterol, mg/dL 


188.4±34.7 


180.7±32.3 


193.6±35.3 


<0.001 


HDL cholesterol, mg/dL 


52.6±13.6 


50.4±14.6 


54.1±12.8 


<0.001 


Triglycerides, mg/dL 


146.1±89.2 


152.5±105.3 


141.7±76.4 


0.114 


25(OH)D, ng/mL 


25.8±12.3 


30.1±13.8 


22.9±10.2 


<0.001 


Glucose, mg/dL 


99 [93 to 108] 


101.0 [94.0 to 112.5] 


97 [91 to 105] 


<0.001 


Insulin, uIU/mL 


7.1 [5.8 to 8.8] 


6.7 [5.4 to 8.5] 


7.3 [5.9 to 9.2] 


<0.001 


HOMA-IR 


1.8 [1.4 to 2.3] 


1.7 [1.3 to 2.3] 


1.8 [1.4 to 2.4] 


0.054 


Hypertension 


295 (36.6) 


128 (39.5) 


167 (34.6) 


0.177 


Diabetes mellitus 


126(15.6) 


65 (20.1) 


61 (12.6) 


0.006 


HOMA-IR >2.5 


170 (21.1) 


60(18.5) 


110(22.8) 


0.172 


Current smoker 


74(9.2) 


68 (21.0) 


6(1.2) 


<0.001 


Current alcohol drinker 


315(39.0) 


192 (59.3) 


123 (25.5) 


<0.001 


Regular exercise 


354 (43.9) 


147 (45.4) 


207 (42.9) 


0.527 



BP, blood pressure; HDL, high-density lipoprotein; HOMA-IR, homeostasis model assessment for insulin resistance. 
Data are expressed as mean±standard deviation, median [inter-quartile range], or number (percent). 
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Table 2. Pearson's Correlation Coefficients of Serum Lipids and 25(0H)D and HOMA-IR (n=807) 







Unadjusted 






Adjusted for sex and age 






HDL 

cholesterol 


Triglycerides 


25(OH)D 


HOMA-IR 


HDL 

cholesterol 


Triglycerides 25(OH)D 


HOMA-IR 


Total cholesterol 


0.354* 


0.223* 


0.006 


0.001 


0.337* 


0.236* 0.063 


-0.012 


HDL cholesterol 




-0.324* 


0.109* 


-0.219* 




-0.320* 0.157* 


-0.231* 


Triglycerides 






0.196* 


0.320* 




0.191* 


0.324* 


25(OH)D 








-0.073* 






-0.054 



HOMA-IR, homeostasis model assessment for insulin resistance; HDL, high-density lipoprotein. 

Triglycerides and HOMA-IR values are log-transformed for analysis. 

*p<0.001. 

V0.01. 

*/k0.05. 



Table 3. Association between Serum 25(0H)D and Log-Transformed HOMA-IR in Total Participants (n=807) 





Unadjusted 


Sex and age-adjusted 


Multiple adjusted 




(3 p value 


P 


p value 


P 


p value 


25(OH)D, ng/mL 


-0.003 0.039 


-0.002 


0.123 


-0.002 


0.247 


Women (vs. men) 




0.056 


0.128 


0.041 


0.272 


Age, yrs 




0.0002 


0.916 


0.001 


0.422 


Body mass index, kg/m 2 








0.065 


O.001 


Current smoker (vs. others) 








-0.026 


0.655 


Current alcohol drinker (vs. others) 








-0.005 


0.877 


Regular exercise (vs. no) 








0.062 


0.052 



HOMA-IR, homeostasis model assessment for insulin resistance. 



sex, age, BMI, smoking status, alcohol intake, regular exer- 
cise, and menopause status, the odds for having insulin resis- 
tance was 1.06 (95% CI, 0.91-1.24). When we examined the 
relation between clinically relevant categories of 25(OH)D 
(<20, 20-30, and >30 ng/mL), we found no association be- 
tween serum 25(OH)D and insulin resistance (Supplemen- 
tary Table 2 and 3, only online). The magnitude and direc- 
tion of these associations were similar to findings in the 
subgroup analyses categorized by smoking status and alco- 
hol intake (data not shown). 



DISCUSSION 



The current study investigated the independent association 
between serum 25(OH)D level and insulin resistance in a 
rural Korean population. Our findings suggest that vitamin 
D is not associated with insulin resistance when adjusting 
for potential confounders in a general population. 

The current literature remains inconclusive with regard to 
whether lower 25(OH)D concentration is independently as- 
sociated with insulin resistance. Several epidemiologic stud- 
ies have reported associations between lower circulating vita- 
min D levels and impaired glucose homeostasis. 1 - 6 " 9 Analysis 



of the U.S. National Health and Nutrition Examination Sur- 
vey datasets revealed that insufficient serum 25(OH)D con- 
centration is significantly associated with type 2 diabetes and 
insulin resistance. 6 As well, lower circulating 25(OH)D lev- 
els were associated with higher fasting insulin and HOMA- 
IR in a middle-aged to elderly Chinese population 9 and a Jap- 
anese working population. 8 The association between lower 
25(OH)D concentration and impaired glucose homeostasis 
was also reported in the Korean population. In the Fourth 
Korea National Health and Nutrition Examination Survey 
(KNHANES), a low serum 25(OH)D concentration was as- 
sociated with diabetes, fasting insulin, and increased HOMA- 
IR. 1 ' 7 However, in the rural Korean population of the pres- 
ent study, we did not observe an independent association 
between serum 25(OH)D level and insulin resistance. Our 
finding is consistent with several previous studies that also 
did not observe an association between serum 25(OH)D 
and insulin resistance. In the Third U.S. National Health 
and Nutrition Examination Survey, an inverse association 
between serum 25(OH)D level and insulin resistance was 
found in only non-Hispanic whites and Mexican Ameri- 
cans, but not in non-Hispanic blacks. 11 In the Hoorn study, 
25(OH)D concentration was not associated with fasting and 
postprandial glucose levels, and incident diabetes. 12 In a 
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healthy Cree community in Quebec, Canada, no association 
between serum 25(OH)D and insulin homeostasis indices 
(HOMA-IR and HOMA-beta) was detected. 14 Moreover, 
even in the Asian studies that reported an inverse associa- 
tion between 25(OH)D and insulin resistance, the associa- 
tion was limited to participants who were overweight or 
obese. 8 - 9 Lastly, a randomized double-blinded trial in wom- 
en with polycystic ovary syndrome and vitamin D deficien- 
cy also did not find vitamin D supplementation to have any 
effect on insulin resistance. 13 

Most of the studies that reported inverse associations be- 
tween serum 25(OH)D level and insulin resistance used na- 
tionwide survey data, and participants were from diverse 
backgrounds. Therein, demonstration of whether the relation- 
ships between serum 25(OH)D level and insulin resistance 
are independent from other risk factors would be difficult, as 
potential confounders would have unlikely been fully con- 
trolled. One of the strengths of our study was that participants 
were relatively homogeneous in ethnicity and geographical 
location. These differences in study populations could be one 
reason for inconsistent findings between our study and the 
nationwide surveys. Another possible reason for this discrep- 
ancy is the difference in the serum 25(OH)D levels. When 
compared with KNHANES results, our study participants 
had relatively high serum 25(OH)D levels; mean serum 
25(OH)D levels were 25.8 ng/mL in our study and 19.7 ng/ 
mL in the KNHANES. Our study was conducted in a rural 
community during the summer, while the KNHANES was 
conducted nationwide throughout the year: serum 25(OH)D 
concentration is known to be significantly higher in rural ar- 
eas than in urban areas and typically higher in the sum- 
mer. 1 ' 711 Vitamin D level is dependent on sunlight exposure 
because ultraviolet light is required for the first step of vita- 
min D synthesis from 7-dehydrocholesterol in the skin. 2 22 

The current study has several limitations. First, this was a 
cross-sectional observation study, therefore a causal relation- 
ship between serum 25(OH)D and insulin resistance cannot 
be made. Second, we did not use the glucose clamp method, 
the gold standard test; instead we used HOMA-IR as a sur- 
rogate measure of insulin resistance. The HOMA-IR equa- 
tion has been widely used in clinical and epidemiological 
studies based on a good correlation between estimates of 
insulin resistance derived from HOMA and from the eugly- 
cemic clamp. 23 25 Furthermore, potential confounding fac- 
tors such as calcium and parathyroid hormone may be re- 
sidual confounders. Third, serum 25(OH)D level was not 
repeatedly measured over four seasons. However, the prima- 



ry aim of the present study was not to estimate absolute lev- 
els of serum 25(OH)D over the year as a descriptive report, 
but to investigate the association between serum 25(OH)D 
and insulin resistance in a general population. If 25(OH)D 
level and insulin resistance are directly related the relation- 
ship should be prominent during the summer season. Thus, 
our findings are not seriously distorted because 25(OH)D 
was measured only during the summer. Last, both serum 
25(OH)D and insulin resistance were measured only once; 
thus, we cannot exclude the possibility that random measure- 
ment error may have weakened the observed association. 

In conclusion, serum 25(OH)D level was significantly as- 
sociated with insulin resistance in the unadjusted model. 
However, the association disappeared after adjusting for po- 
tential confounders, such as sex, age, BMI, smoking status, 
alcohol intake, and regular exercise. This study suggests that 
vitamin D is not independently associated with insulin resis- 
tance in a rural-dwelling Korean adult population. 
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